Mesoscopic conductance fluctuations in YBa 2 Cu 3 07_j grain boundary Junction at low 

temperature 
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The magneto-conductance in YBCO grain boundary Josephson junctions, displays fluctuations 
at low temperatures of mesoscopic origin. The morphology of the junction suggests that transport 
occurs in narrow channels across the grain boundary line, with a large Thouless energy. Nevertheless 
the measured fluctuation amplitude decreases quite slowly when increasing the voltage up to values 
about twenty times the Thouless energy, of the order of the nominal superconducting gap. Our 
findings show the coexistence of supercurrent and quasiparticle current in the junction conduction 
even at high nonequilibrium conditions. Model calculations confirm the reduced role of quasiparticle 
relaxation at temperatures up to 3 Kelvin. 



I. INTRODUCTION 

In high critical Temperature Superconductor (HTS) 
junctions, including grain boundary (GB) structures, 
it is well established that various interplaying mecha- 
nisms contribute to transport with different weights still 
to be completely defined in a general and consistent 
framework 1 ' 2,3 . As in all non-homogeneous systems, the 
barrier region will significantly contribute to determine 
the transport properties across the structure. What is pe- 
culiar of HTS is the complicate material science entering 
in the formation of the physical barrier microstructure. 
This will depend on the type of device and the fabrica- 
tion procedure. The material science complexity of HTS 
may also result in different precipitates and inclusions 
present at interfaces and grain boundaries, and in some 
type of inherent lack of uniformity of the barriers 2 -. All 
this has turned into some uncertainty about the nature 
of the barrier and has originated various hypotheses on 
the transport properties. The most wide-spread models 
are basically all in-between two extreme ideas 2 -^^^: 
on the one hand resonant tunnelling through some kind 
of dielectric barrier— on the other, especially in GB 
junctions, a barrier composed of thick insulating regions 
separated by conducting channels, which act as shorts 
or microbridges 8 -^. In most cases the interface can be 
modeled as an intermediate situation between the two 
limits mentioned above. A transition from one extreme 
to the other can therefore take place. What is un- 
fortunately missing is a way to describe this tuning- 
transition through reliable and well defined barrier pa- 
rameters (for instance the barrier transparency). The 
predominant d-wave order parameter symmetry (OPS) 



is another important 1 , to which a large part of the phe- 
nomenology has been clearly associateoU* 2 -^ ' 10 ' 11 . D-wave 
OPS implies the presence of antinodal (high energy), 
and nodal (low energy) quasi-particles in the conduction 
across junctions and the absence of sharp gap features in 
the density of states of the weak link. Recently, low tem- 
peratures measurements have proved macroscopic quan- 
tum tunneling (MQT) in YBCO GB junctions, stimulat- 
ing novel research on coherence and dissipation in such 
complex system o 12 ' 13 ' 14 . 

In this work we report on an investigation of magne- 
toconductance at low temperatures for the same type of 
biepitaxial GB junctions^ 5 ., used for the MQT experi- 
ments. These structures are very flexible and versatile, 
guaranteeing on the one hand low dissipatio n 13 ' 14 and 
on the other a reliable way to pass from tunnel-like to 
diffusive transport on the same chip by changing the in- 
terface orientatio n 12 ! 15 . We give direct evidence of the 
role played by narrow conduction channels across the GB. 
These channels may have different sizes and distributions 
and obviously a different impact on the transport prop- 
erties. When increasing applied voltage, mesoscopic con- 
ductance fluctuation s 16 ' 17 ! 18 ' 19 ! 20 appear in our samples, 
at low temperatures, not dissimilar from what is usually 
observed in normal narrow metal samples 2 —. We expect 
that, in our sample, typical sizes of the current-carrying 
constrictions L x and L y range from 50 nm to 100 nm 
and, as a consequence, Thouless energy E c ( see Table 
U ) turns out to be quite large, when compared to the 
values usually measured in traditional normal metal ar- 
tificial systems^. The mesoscopic effects persist at volt- 
ages about twenty times larger than the Thouless energy. 
" Novel" mesoscopic issues that emerge from the analyses 
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carried out in the present work are tightly connected to 
the nature of the GB systems: a) a smooth crossover ap- 
pears to exist from the coherent conduction mostly driven 
by the supercurrent, to the magnetoconductance driven 
by quantum coherent diffusion of quasiparticles across 
the mesoscopic area, when the voltage at the junction in- 
creases; b) in analogy to pairs, quasiparticles also appear 
to have a large phase coherence as proved by the shape 
of the power spectrum of the conductance fluctuations, 
up to temperatures of 3K; c) the voltage drop appears 
to be concentrated at the GB, and non equilibrium does 
not affect substantially the mesoscopic interference over 
a wide area about 1/im 2 . 

This work builds upon a previous report, where the 
main ideas have been illustrated 23 . Herein a more com- 
plete analysis of the experimental data is carried out. 
We have applied the " protocol" established in the last 20 
years on semiconducting and normal metal nanostruc- 
tures to our system and we have extracted the character- 
istic lengths and scaling energies. 

In Section II we give some details about the sample 
fabrication. By presenting the magnetic mesoscopic fin- 



gerprints of our sample in Section IIIA, we collect ev- 
idences of the mesoscopic character of the conductance 
fluctuations that we have measured. In section IIIB we 
derive from the ensamble average of the fluctuations the 
variance of the conductance, which is presented in section 
IIIC. 

In Section IV we show the conductance autocorrelation 
for different magnetic fields in an intermediate voltage 
range. By analyzing the power spectral density we esti- 
mate the phase coherence length L v which is found to be 
< L, at intermediate voltages ( V — 14 ~ 18 mV). The 
conductance autocorrelation at different voltages allows 
us to interpret the role of noncquilibrium by defining the 
voltage dependence of the phase coherence length, as dis- 
cussed in Section V. The discussion of the results can be 
found in Section VI. Our simple model theory well ac- 
counts for the experimental results and clarifies the sur- 
vival of non locality in the quantum diffusion in presence 
of a large voltage bias. Table U gives additional informa- 
tion on the planning of this work. Our conclusions can 
be found in Section VII. 



TABLE I: Summary of the results: 



Topic 


Relative measurement 


Extracted parameters 


Estimates 


Junction geometry 


Magnetic pattern: Ic vs. H in Fig.Q 


L x ~ L y ~ 50nm 


E c ~ ImV 


Transport properties 


Rn in Fig.gJ) and I C R N in Fig.© 


Rn = 400fl, T > 257mK 


I C R N ~ E c /e 


Mesoscopic Fingerprints 


Resistance fluctuations in Fig.s(l5l6l) 


var[g] < 1 




Phase coherence length 


Autocorrelation vs AH in Fig.Q, PSD in Fig.sjHEJl 






Coherent phase breaking time 


Autocorrelation vs AV in Fig. (1101) 


t v ~ 400 ps 





II. FABRICATION AND AVERAGE 
TRANSPORT PROPERTIES OF THE SAMPLE 



The GB Josephson junctions are obtained at the inter- 
face between a (103)YBa 2 Cu 3 7 _ s (YBCO) film grown 
on a (110) SrTiOs substrate and a c-axis film deposited 
on a (110) Ce02 seed layer (see Fig.([T])). The presence 
of the Ce02 produces an additional 45° in-plane rotation 
of the YBCO axes with respect to the in-plane direc- 
tions of the substrate 1 ^. The angle 8 of the grain bound- 
ary relative to the substrate a, b axes is defined by suit- 
ably patterning lithographically the Ce02 seed layer (see 
Fig. ((T]a)). Details about the fabrication process and a 
wide characterization of superconducting properties can 
be found elsewher e - 12 ! 15 . The interface orientation can be 
tuned to some appropriate transport regime, evaluated 
through the normal state resistance Rn and critical cur- 
rent density (3c)- Typical values are reported in Fig.@ 



and compared with data available in literature 3 -. 

In the tilt cases 3c ~ I0 3 A/cm 2 and gn—^KRnA) 
fa 0.2(m^cm 2 )" 1 , both measured at T = 4.2 K (A is the 
junction cross section). Twist GBs junctions are typically 
characterized by higher values of 3 c in the range 0.1-4.0 
x 10 5 A/cm 2 and a N »1Q (mficm 2 )" 1 (at T = 4.2 K). 

We have selected YBCO grain boundary junctions and 
measured their I-V curves at low temperatures T, as a 
function of the magnetic field H applied in the direction 
orthogonal to the plane of the junction. In HTS junc- 
tions, the correlation between the magnetic pattern and 
the current distribution profile along the junction is made 
more complicate by the d-wave order parameter symme- 
try (OPS), which generates anomalous magnetic response 
especially for faceted interfacesi 2 -^^. Additional devia- 
tions are expected because of the presence of the second 
harmonic in the current-phase relatio n 10 ' 11 . In Fig.© 
we report the magnetic field dependence of the IV char- 
acteristics of the maximum Josephson current of the of 
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FIG. 1: a) Sketch of the YBCO biepitaxial GB junction used 
in this experiment. b)Geometry of the model system of the 
current-carrying constriction. 



FIG. 3: (color online) The I-V characteristics as a function of 
magnetic field. 




FIG. 2: (color online) The product I c Rn is reported as a 
function of the critical current density J c . Data are collected 
from Raf.8LUUU8.ai 



the junction that we have extensively investigated in this 
work (with barrier orientation 0=60°). This angle gives 
the maximum Jo^- The magnetic response presents a 
maximum of the critical current at zero field and two al- 
most symmetric lobes for negative and positive magnetic 
fields respectively. At higher magnetic fields (above 100 
G or below -100 G), the critical current is negligible. The 
flux periodicity is roughly consistent with the size we ex- 
pect for our microbridge (50-100 nm), since the London 
penetration depth in the off-axis electrode is larger than 
the one in c-axis YBCO films, of the order of microns 
(see Refjil and^ for instance). Experimental data can 
be compared with the ideal Fraunhofer case in the crud- 
est approximation without taking into account the pres- 
ence of a second harmonic or any specific feature of HTS. 
Even if deviations from the ideal Fraunhofer pattern are 
present, they can be considered to some extent minor if 
compared with most of the data on HTS grain boundary 
Josephson junctions, which present radical differences. 
We can infer an uniformity of the junction properties 



approximately on an average scale of 20-30 nm, which 
is remarkable if compared with most results available in 
literature. Even if we cannot draw any conclusion on 
the current distribution on lower length scales, we can 
rule out the presence of impurities of large size along the 
width of the active microbridge. In fact, were there more 
than one active microbridge, the current of each of them 
would add in parallel and the pattern would present other 
periodicities referring to the area enclosed between the 
conduction channels (Ref. 30 ). The I-V characteristics of 
the HTS Josephson junctions still present features, which 
cannot be completely understood in terms of the classi- 
cal approaches used to describe the low critical temper- 
ature superconductors Josephson junctions. These are 
frequently observed and often referred in literature as un- 
conventional featuresi^. Examples are 3 -: a) IcR-w val- 
ues are much lower than the gap value A; b) the shape of 
the I-V strongly depends on the critical current density; 
c) I-V curves show significant deviations from the Resis- 
tively Shunted Junction model (RSJ); d) there is a poor 
consistency between the amplitude of the hysteresis and 
the extracted values of the capacitance, when compared 
to low-Tc superconductor junctions. 

At low temperatures, the resistance vs applied current 
R(I), as derived from the I-V characteristics is rather 
temperature insensitive, while the critical current Ic 
maintains a sizable temperature dependence. In the inset 
of Fig. (J4j) we show the resistance R(I) at zero magnetic 
field for three temperatures: 257mK, IK and 3K. R(I) 
obviously vanishes in the Josephson branch and displays 
a sharp peak when switching to/from the finite voltage 
conductance. The data are displayed in order to show 
the critical current at I > and the retrapping current 
at I < 0. In Fig. (|4]main-panel) we show a blow up of 
the resistance R(V) in a range of voltage values V be- 
tween 0.25 mV and 30 mV at zero magnetic field and for 
different temperatures T = 257mK, IK, 3K. 
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Measurements have been taken after different cool- 
downs in the time lapse of two years to study the sample 
dependent properties. 

The average resistance in the range of voltages V ~ 
IQmV ^rlhmV has been stable for about eighteen months 
at ~ 180 il and has increased in the last year up to about 
430fi. These changes should be attributed to aging of the 
diffusion properties at the grain boundary. However, in 
the meantime, no significant change of the Ic has been 
detected. Only one sample was available with such a 
reduced width. The steady progress in nanotcchnology 
will probably lead to the realization of reliable micro- 
bridges of nominal width of a few hundred nanometers, 
from which it will be easier to have junctions with trans- 
port carried by very few mesoscopic channels. We finally 
signal a strong similarity of the I-V and dl/dV-V curves 
of our junctions with those from sub-micron YBaCuO 
junctions, reported in Refl3ll 
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FIG. 4: (color online) Main: zooming on the oscillations 
of R vs. I above I w 80fiA. Inset: Resistance as a func- 
tion of applied bias current for three different temperatures: 
257mK, IK, 3K. 



III. MAGNETORESISTANCE AND 
MESOSCOPIC FINGERPRINTS 

A. Resistance fluctuations 

According to what reported in the previous Section, 
we figure out that most of the current in the junction 
substantially flows across a single nanobridge of charac- 
teristic size L x < lOOnm. Possible lack of spatial uni- 
formity of the current distribution on a scale of less than 
20 nm does not affect the arguments developed below. 

We call x the flow direction and y the direction per- 
pendicular to the nanobridge, as shown in Fig.fljD). 
The diffusion coefficient D in YBCO is expected to be 
~ 20 -7- 24 cm 2 / se<f£. By estimating the Fermi velocity 



vp w 7x 10 7 cm/ sec for optimally doped YBCO, we con- 
clude that the mean free path £ is smaller than the size 
of the nanobridge. Therefore we argue that the trans- 
port in the junction is diffusive, which is confirmed by 
the observation of the resistance fluctuations. 

The Thouless energy, as derived from the expected size 
of the nanoconstriction, is E c = hD/L 2 > 1 meV. This 
value is confirmed by our measurements as discussed in 
Section III. The number of transverse scattering channels 
in the constriction for a fixed cross section A is approx- 
imately M ch w k 2 F A ~ 5 ■ 10 4 . A ~ 100 x 100 nm 2 is 
given by the product of the thickness of the film and 
the width of the channel. Hence, quantization of trans- 
verse levels in the bridge does not seem to play any role 
even at the lowest temperatures investigated. Indeed, 
ksT » 5 ~ E c /Mch ~ O.lfieV, where 5 is the mean 
energy level spacing. As a consequence, our system can 
be thought as a disordered bridge in the diffusive limit 
i < L x . 

We concentrate on the marked non periodic fluctua- 
tions of the resistance at finite voltages, with magnetic 
field in the range B = —100 -j- 100G. An example of the 
magnetoresistance fluctuations is reported in Fig. ([5]) top 
panel for T = 257mK, IK, 3K. 

The fluctuations are not related to the magnetic depen- 
dence of the Josephson critical current I C (H) at voltages 
V > 2mV. 

In order to avoid trapping of flux which may oc- 
cur when increasing the temperature, especially at the 
higher fields, the data shown here refer to a single cool- 
down. Occasionally the pattern has still a slight de- 
viation from reproducibility within one single cooling 
bath, which could be due to finite relaxation in the 
spin orientation of paramagnetic impurities. CONFER- 
MARE The resistance pattern derived from our four ter- 
minal measurement does not show any mirror symmetry 
Rn(H) ^ Rn(—H). Below IK there is little temperature 
dependence. The amplitude of the fluctuations decreases 
between T = IK and T = 3K. They are sample depen- 
dent, as different cool-downs provide different patterns. 
All these features, as well as the ones described below, 
strengthen the conclusion that they are mesoscopic fluc- 
tuations. 

The color plot of the resistance fluctuations in the I, H 
plane provides the fingerprints of our sample. The devi- 
ation from the average SR = R — R is shown in Fig. ([5|) 
(three color-plot panels) as a color plot (color online), 
for three different temperatures, 257 mK, IK, 3K. R is 
the average resistance performed over the full range of 
magnetic fields. The pattern keeps its shape within one 
single cool down and the contrast of the colors increases 
in lowering the temperature. The color scale is such that 
the dark red color refers to resistances significantly larger 
than the average, while the dark blue color refers to resis- 
tances significantly smaller then the average. The data 
have been filtered by gaussian convolution to get rid of 
the underlying white noise. 

Despite the small equilibrium thermal length Lt — 
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FIG. 5: (color online) Color plot of the resistance fluctua- 
tions 8R(H, I) as a function of the applied current and of the 
magnetic field for three temperatures: 257mK, IK, 3K. The 
right panel shows SR(H — 2.7G, /) vs. I for the three differ- 
ent temperatures. The three top panels show single magneto- 
conductance trace for each temperature at fixed bias current 
(indicated in the labels). 

yjTiD/k B T ~ 0.14 turn at 1 K, there is a clear persistence 
of the fingerprints up to T = 3K. This suggests that the 
strong non-equilibrium conditions induced by the applied 



voltage do not allow the thermalization of the carriers in 
the sample. Indeed, our results do not change qualita- 
tively up to 1.5 K , and transport can be classified as non 
equilibrium quantum diffusive because L v < L x . Here 
L v is the phase coherence length for carriers diffusing in 
the junction area and is < lfim (see Section IV). 



B. Ensemble average 

The variance war[g] is the ensemble average of the am- 
plitude squared of the conductance fluctuations, ((Sg) 2 ), 
where g = G/(2e 2 /h) is the dimensionless conductance. 
We analyze the fluctuations of the conductance obtained 
by averaging over runs at different magnetic fields up to 
100G, at different voltages. Here we argue that this av- 
erage can be taken as an acceptable ensemble average 
and provides a bona fide information about the variance 
i>ar[g] of the conductance and its autocorrelation. To jus- 
tify our statement, we have to show that the Cooperon 
contribution to the variance is not significantly influenced 
by H up to at least 100 G. The variance of the conduc- 
tance at equilibrium at temperature T can be calculated 



var[g(H,T)} 
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where s is the spin degeneracy, f(x) = (xcothx — l)/sinh 2 x and Fd,c are the Diffuson and Cooperon autocorrelation 
functions, respectively. They can be rewritten in terms of the eigenvalues A„ ,c of the diffusion equation^!: 
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Here Dt = P/d ~ 4 • 10 cm (where d is the effective dimensionality ) and Ti„ is the inelastic relaxation time 
(nn » t ~ 0.2 psec). Ai,A 2 are the vector potentials (Hi and H 2 the magnetic fields) influencing the outer and 
inner conductance loops respectively and the +/— sign refers to the C/D propagator respectively. We have: 



F D , C = L^(Dr) 2 £ 



1 



X D,C 



5Re- 



1 



Act 



(3) 



At zero temperature only AE = contributes to the integral in Eq.([T]), so that the eigenvalues become real. In the 
evaluation of the variance, Hi — H 2 — H implies that the Diffuson eigenvalues become insensitive of the magnetic 
field. Instead, the Cooperon eigenvalues depend on 2H and can be written in analogy with the Landau levels energies. 
It follows that : 
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where u>h = AeHD /Tic. Eq.Q defines a decay threshold from the truncation of the sum over the orbital quan- 
field of the Cooperon H Q ~ irhc/ (\2eDr), which derives 
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turn number n at n max as %/mcj c l 2 (uj c = eH/mc is 
the cyclotron frequency). This limitation is required by 
quantum diffusion ((r 2 )„ max > I 2 ). The large value of 
Timaxf In our case (~ 10 2 for H < 600 G), determines 
H a > 2.5 Tesla which is far beyond the field strengths 
that can be applied to our sample without trapping flux 
due to vortices. This confirms that averaging over the 
interval of H values H G (—100 -r- 100 G) is equivalent 
to a sample average, without introducing significant field 
dependencies. In the following the ensamble average will 
be denoted by the symbol (...}#. As it is shown in the 
next Section, the typical magnetic field scale that arises 
from the autocorrelation of the conductance is ~ 10 G, 
much smaller than the interval over which the average is 
performed. 

In the rest of the paper, we will generically denote 
the conductance autocorrelation, which is an extension 
of Eq.([I]). by K g . This quantity depends on many vari- 
ables: T,H = (Hi + H 2 )/2,AH = H x - H 2 ,V = 
(Vi + V 2 )/2,AV = Vi — V 2 . When no ambiguity arises, 
we have taken the liberty to list just the parameters rel- 
evant to the ongoing discussion, in order to simplify the 
notation. 



C. Variance of the conductance and different 
voltage regimes 

The conductance is derived from the I/V character- 
istic. We have checked the behavior of the differential 
resistance, measured through a standard lock-in method, 
and we have found qualitatively similar results. 

In Fig.® (upper panel ) we have reported the mea- 
sured conductance fluctuations vs voltage bias and mag- 
netic field at the temperature T = 257mK, in a grey scale 
plot. The plot shows two different regimes: 

a) Low voltages ( V < 3mV), where fluctuations ap- 
pear to be very high. Fluctuations in this range 
mostly arise from precursive switching of the cur- 
rent out of the zero voltage Josephson state. The 
analysis of this range of voltages is better dis- 
cussed within the Macroscopic Quantum Tunneling 
dynamics^; it requires full account of the super- 
conductive correlations and is not addressed in this 
paper. 

b) Large voltages (V > 5mV) . In this regime we ob- 
serve some reproducible, non periodic and sample 
dependent fluctuations. The variance of the con- 
ductance < (Sg) 2 >h, is plotted vs voltage bias ( 
Fig. ([6]) bottom panel) for two temperatures. The 
scale for its magnitude is estimated according to 
5R/R = Sg/g w 0.01 with R = 410f2(T = 300mAT). 
The variance v ar[g] stabilizes around unity at V > 
7mV. As the voltage increases V > 18 mV, the 
variance is increasingly reduced. However, small 
amplitude fluctuations seem to persist over a wide 
voltage range up to values which are by far larger 



than those in normal constrictions. Fluctuations 
survive up to voltages which are many times the 
Thouless energy. 

Here we focus on the variance var[<?] and on the auto- 
correlation of the magnetoconductance as a function of 
voltage at low temperatures up to ^ >> E c /e ~ 1 mV . 
Our data can be interpreted on the basis of models for 
the quantum interference of carriers transported in the 
narrow diffusive channel across the GB line. The ex- 
perimental findings are consistent with a large Thouless 
energy E c and quite long dephasing times t v . A com- 
parison of the data with the results of our models seem 
to confirm that non equilibrium effects induced by the 
voltage bias V are not the source of heavy energy relax- 
ation of the carriers, even at voltages V >> E c /e. A 
discussion about the voltage dependence of the variance 
of the conductance for large voltages can be found in the 
Section VI. In the next Section we analyze the conduc- 
tance autocorrelation at finite voltage in some detail, to 
extract information about the phase coherence length L v 
and the phase coherence breaking time t v . 

IV. SAMPLING NONLOCALITY: 
AUTOCORRELATION vs AH 

The variance var[<?], of the conductance fluctuations 
((8g) 2 )H discussed in Section III.C can be derived from 
the maximum at AH « of the more general autocorre- 
lation function: 

K g (V, AH) = (5g(V, H + AH)Sg(V, H)) H (5) 

The data have been averaged over H, as usual, as well as 
over a small interval of voltage values about V . K g is the 
sum of Cooperon Fc and Diffuson Fry contributions. Due 
to the independence of our results of the temperature we 
consider the zero temperature limit of Eq. ([T]) . In Fig. ([7]) 
we plot the measured autocorrelation of the dimension- 
less conductance vs. AH at T = 257mA" for three values 
of V: 10 mV (blue curve), 14 mV (red curve) and 18 mV 
(red curve) (color online ). The data have been averaged 
over a voltage interval of width SV = 0.5mV (the results 
do not depend on this choice). 

Similar curves have been measured for normal metal 
wires at zero voltage bias 2 ^. 

The autocorrelation of Fig. (JTJ) is practically insensi- 
tive to increasing temperature up to about 1.5 K. This 
fact can be viewed as evidence that a significant contri- 
bution to transport and to the conductance fluctuations 
is still provided by the pair current. Its time average 
and absolute value can be seen as rather temperature 
and voltage independent, at fractions of Kelvin, while 
quasiparticles remain rather frozen, provided the voltage 
does not increase too much—. For weak links, charac- 
terized by and higher barrier transparency the contribu- 
tion of the contribution of the supercurrent in the I-V 
curve can be relevant at finite voltages^. The physical 
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FIG. 7: (coJor online) Autocorrelation < (Sg) 2 > vs AH at 
T = 257mA' for three different values of the voltage drop: 
V = 10 (black curve), V = 14 mV (red curvej, V = 
18 mV" (green curve). Inset. Results of the model calculation 
sketched in the text for £ = 1.05 (black curve), 1.15 (red 
curve), 1.4 (green curve). 
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FIG. 6: (color online) Top panel: grey-color plot of the fluctu- 
ations of the dimensionless conductance as a function of the 
voltage V and of the applied magnetic field H at T = 257mK. 
Bottom panel: Variance of the dimensionless conductance as 
a function of the voltage V around zero magnetic field for 
T=257mK and T=3K. 



reason is that the phase changes in a sharply nonlinear 
manner with the greater part of the period being close to 
7r/2 + 2wr. In addition, non-equilibrium effects 3 ^ and un- 
conventional order parameter symmetry (with a not neg- 
ligible second harmonic component in the current-phase 
Josephson relation^) are possible additional sources of 
supercurrent flowing at finite voltage. Our results seem 
to confirm the presence of non negligible contribution of 
supercurrent at large voltages, from a different perspec- 
tive. This is consistent with Ref.s [36ll37ll38l and possibly 
with the observation of fractional Shapiro steps on YBCO 
grain boundary Josephson junction a 39 i 40 . 

According to the remark made above, we can assume 
that the current at T ~ is only a function of the phase 
difference ip between the two superconducting contacts. 
This assumes little dephasing induced by inelastic scat- 
tering processes, but not necessarily the absence of quasi- 



particle contributions to the current which still depends 
on ip. 

The inset of Fig. ([7]) shows the result of a simple model 
calculation of the autocorrelation function based on the 
following assumptions: 1) negligible proximity effect in 
the sub-micron bridge induced by the superconducting 
contacts; 2) an equilibrium approach to transport, in 
which the current is mostly phase dependent; 3) handling 
of the magnetic field H is treated as a small correction 
and therefore H only appears in the gauge invariant form 
of the phase difference. 

The model (see Ref4i for details) uses as unique fit- 
ting parameter £ = L y /L v , where L y is the transverse 
size of the conduction channel. The curves plotted in the 
inset are with £ = L y /L v =1.05 (black curve), 1.15 (red 
curve), 1.4 (green curve). The three different measured 
curves in the main plot of Fig. (J7|) refer to different bias 
voltages and cannot be directly compared to the theo- 
retical curves in the inset of Fig.(J7J). The qualitative 
agreement between experimental and theoretical curves 
is evident, provided we assume that £ increases, with in- 
creasing voltage. This assumption is feasible since, on 
the one hand L v is likely to be reduced when increasing 
applied voltage, and the number of conduction channels 
increases by changing the voltage and, as a consequence, 
the effective width of the bridge. If we assume that L v 
scales with voltage as V^ 1 ^ 4 ( see discussion in Section 
V ), we find that the £'s, that have been chosen to draw 
the inset of Fig. (J7J), are consistent with the voltages of 
the experimental curves within 15% of error. 

The qualitative fit, based on the simple theoretical 
model used here, gives evidence of the fact that non- 
equilibrium does not seem to spoil the autocorrelation as 
a function of the magnetic field, even if the voltage bias 
exceeds the Thouless energy: eV > E c . While the de- 
phasing time t v is discussed in the next Section, here we 
are in position to extract the value of the phase coher- 
ence length L v from the Power Spectral Density (PSD) 
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of the conductance autocorrelation function. 




f (gauss 1 ) 



FIG. 8: (color online) PSD at T= 257 mK, IK, 3K and V = 
7m V , averaged over a voltage interval of 0.5mV . The curves 
have been shifted for clarity. 

The Power spectral density (PSD) of the conductance 
autocorrelation is plotted in Fig.© vs fn, the conjugate 
variable to the magnetic held AH , for data at T = 273 
mK, IK, 3K and V= 7 mV. Jh has the dimension of an 
inverse magnetic held. Inspection of Fig.© shows that 
there is a linear slope at small frequencies and a roughly 
flat trend at larger frequencies. The latter is due to the 
white noise affecting the measurement. Curves have been 
shifted, in the hgure, for clarity. In reality, both the linear 
slope and the constant value are almost independent of 
the temperature. The linear slope allows to extract the 
value of the phase coherence length L v according to the 

log(PSD/PSD(f H = 0)) = -2nf H H c + cnst , (6) 

as correlation held H c can be related to L v as follows: 
L v ~ y / hc/2eH c . From this plot we derive H c ~ 10 G, 
that gives L v < f fim. A logarithmic plot of the V de- 
pendence of L V (V) is reported in fig©. The straight line 
drawn among the experimental points shows the func- 
tional dependence V -1 / 4 . 



V. NON EQUILIBRIUM EFFECTS ON 

In this Section we study the dependence of the auto- 
correlation function on V and AV: 

K g (V, AV) = (Sg(V + AV, H)6g(V, H)) H . (7) 

A voltage average has been performed over inter- 
vals centered at V, of typical size lower than the 
min{E c , AV}. 

In Fig. (|1 Oh ) we report the experimental autocorrela- 
tion function K g (V, AV) vs AV for various values of the 



FIG. 9: Logarithmic plot of L v vs V derived by using the 
eq.©. Full lines are a guide to the eye. The slope of the 
straight line is —1/4. 
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FIG. 10: ( coJor online) The measured autocorrelation of 
the conductance K g (V, H — 0, AV) vs. the voltage difference 
AV, for average voltages V ~ 7.5, 12.5, 17.5, 25.0 mV ( panel 
a)); theoretical fit based on Eq. (|10[) (panel b)); Resistance 
R vs. V (panel c)). There is correspondence between the 
colored ranges of panel c) and the colors of the curves of the 
panels a) and b). 



applied voltage V . We stress that the applied voltage is 
larger than the two natural energy scales: E c /e ~ 1 mV 
and the nominal superconducting gap A/e ~ 20 mV. 
The tail of the curves shows a large anti-correlation dip 
when AV increases and damped oscillations. The auto- 
correlation maximum flattens when the applied voltage 
V increases. Fig.([TU]c) emphasizes the voltage range to 
which each curve of the panel a refers, by using the same 
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color. 

We have reproduced the same trend of the data in 
FigJlO] a), by assuming that proximity effect induced by 
the superconducting contacts does not play an important 
role in the transport and by adapting th e quas i- onc- 
dimensional non-equilibrium theory of Ref. (|43l)44T ) to our 
case (see Fig fTUl b)). We will report on the derivation of 
our theoretical results elsewhere^. We just mention that, 
our nonequilibrium approach gives rise to an autocorre- 
lation K g (V, AV) which is a function of x = eAVL 2 /hD 



and of the parameters t v /t c /d an d (L/L v ) 3 . Here 
r c j D = e 2 L 2 y D{H l ± H 2 ) 2 /{\2h 2 c 2 ) (with the + (-) sign 
for the C (D) case ) is the relaxation time induced by 
the magnetic field^. In the limit of vanishing AV, the 
result for the variance is recovered, which, up to numer- 
ical factors, is given in terms of the ratio between the 
Airy function Ai(u) and its derivative with respect to 
the argument Ai'(u): 



t v ,L v and t v /t u all depend on the voltage. This for- 
mula is similar to the thermal equilibrium result by Alt- 
shuler, Aronov, Khmielnitski (AAK) 46 , which includes 
the dephasing induced by e-e scattering with small en- 
ergy transfer. In particular AAK find 



1 



7~(p Lip 



e 2 k B T 
h 2 a 



(9) 



where a is the conductivity. At intermediate applied volt- 
ages, our zero temperature prefactor Dt v L v /L 3 in Eq.© 
looks similar to the AAK prefactor TiDL v /(3-KL 3 kBT) if 
Ti/t v replaces fcsT. 

By choosing an appropriate voltage dependence of the 
fitting parameter t v /t u and L v /L, we find that the au- 
tocorrelation scales with L V (V) as follows: 



K g {V,AV)^ T 



v=C.D 



AV 

V c ' 



L 



(10) 



with Dt, 



T J 2 



and Tc — tq, V c = h/r v . The explicit 
form of the function T will be given elsewhere^. Its 
limiting form for AV — ► gives Eq.©. This scaling law 
is exploited to plot the curves of Fig pOb V 

The scaling among the blue, red and cyan curves re- 
produces reasonably well the experimental pattern. This 
indicates that the V dependence of eq. fTt)]) is well ac- 



counted for by simply reducing L v with increasing V. 
On the contrary, the black curve at the lowest voltage 
V ~ 7.5mV requires adjusting the prefactor after the 
scaling, to make the central peak higher and narrower. 
This could be a hint to the fact that, when the voltage is 
rather low, the superconducting correlations may be rel- 
evant and should be included in deriving the functional 
form of the T function. Needless to say, an increase of 
the t v /t v parameter implies a reduction of the value of 
var[g). 

In writing Eq. (|10|) . extra contributions arising from 
the non linear response have been neglected. Actually, 



Eq. (|10p would be the result within linear response theory 
only, except for the V dependence of L V (V). 

The global interpretation of the data given here shows 
that a monotonous decrease of r v with increasing voltage 
is not achieved. A non monotonous decrease of var[g] vs 
V is indeed found, as can be seen from Fig.©. According 
to the correspondence Dt v = L 2 , and to Fig.© a general 
decreasing trend of K g (V, AV) with increasing voltage 
could take over only above 20 mV. In Eq.© derived from 
our model calculation, we have found h/r v in place of 
fcgT appearing in the AAK result of Eq.©. At larger 
voltages the expected substitution in Eq.© is^ ksT — > 
eVL v / ' L and, by requiring the consistency, 



L 



t P\eVL v /L 



L 



D_ hg_ 

L 2 eV 



(11) 



where g = ha/(e 2 L) in 1 — d. This would give a decay 
law for the coherence length L v ~ V^ 1 ^ 4 , which is not 
clearly recognizable in our experiment. 



VI. DISCUSSION 

We recollect here the main experimental facts that can 
be extracted from our data regarding quantum transport 
in a GB YBCO JJ . 

The magnetic dependence of the maximum critical cur- 
rent suggests an active transport channel of the order 
of 50-100 nm. Uniformity of the critical current is on 
scales larger than about 20 nm. The Thouless energy 
E c = TiD/L 2 turns out to be the relevant energy scale in 
this case. The normal resistance Rn of the HTS junction 
is of the order of 200f2, increasing up to 480f2 with time, 
due to aging of the sample. The zero field Josephson 
critical current appears to satisfy IqRn ~ E c /e. This 
product is definitely much smaller than A/e, where the 
nominal superconducting gap is A = 20 meV. This rep- 
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resents additional evidence that the proximity effect in- 
duced in the bridge in the absence of applied voltage is of 
mesoscopic origin. The superconductive pair coherence 
length £ s < L v < L, as opposed to the classical regime, 
L < £ s , when the tail of the order parameter enters both 
superconductors of the junction 4 ^. We speculate that the 
oscillations in the resistance as a function of V, shown in 
the inset of Fig. ([!]), could be due to this mesoscopic ori- 
gin. 

Remarkable conductance fluctuations have been found 
in a voltage range up to 20E C in the magnetic field range 
of H e (-100, 100) G for temperatures below 3 K. In the 
explored window, we do not measure an halving of the 
variance with increasing field H^. The crossover field 
H Q at which the Cooperon contribution to the variance 
is expected to disappear, as given by Eq.((3|) is estimated 
of the order of few Teslas 4 ^. Decoherence induced by the 
Zeeman energy splitting requires even larger fields. 

Transport has been measured in highly non- 
equilibrium conditions. Hence the temperature depen- 
dence is quite weak up to T ~ 1.5 K. We have concen- 
trated our analysis in the voltage range V £ (7, 30) mV, 
where the conductance fluctuations reach a steady value 
for the variance at T below 1 Kelvin, var[g](V) < 1. 
These properties confirm that the fluctuations are due to 
quantum coherence at a mesoscopic scale. 

The PSD of the autocorrelation of the conductance 
at different fields AH allows to identify H c as 10 G as 
the field scale for the mesoscopic correlations, weakly 
dependent on V. This value leads to a phase coher- 
ent length L v < 1 fxm. We have plotted the coher- 
ence length extracted from the autocorrelation PSD vs 
V in Fig!] and compared it with L V {V) ~ (V/V )- s 
with s ~ 0.25^. A similar exponent has been found 
in a limited range of voltage bias V > E c /e in gold 
samples 49 , in which UCF (universal conductance fluctu- 
ations) and Aharonov-Bohm oscillations were found. De- 
phasing mechanisms are low frequency electron-electron 
interaction, magnetic impurity-mediated interaction^ 
and non equilibrium quasiparticle distribution 51 . Ac- 
cording to FigEl the comparison is not conclusive. As 
a matter of fact all data of conductance autocorrelation 
at finite voltage reported for normal wire s 22 ! 49 identify a 
Thouless energy E c ~ 1/ieV, three orders of magnitude 
smaller than in our HTS device and refer to applied volt- 
ages not larger than ml^'s. Still, mesoscopic coherence 
persists in our sample, up to voltages much larger than 
the Thouless energy. We do not find any linear increase 
of the conductance autocorrelation with voltage at large 
voltages^ 3 -. 

Our model calculation appears to reproduce the gross 
features in the dependence of the conductance autocor- 
relation K g (V, AH, AV) on AH as well as on AV. In 
the case of K g (V, AH, 0) we limit ourselves to the linear 
response term only and the effect of the voltage bias just 
appeared as a small reduction of L V (V) with increasing 
V. 

To model the trend of K g (V, 0, AV) vs AV given 



by the experiment, non-equilibrium cannot be i gno red. 
Our derivation extends the calculation of Ref.|43.44). 
We give a simple estimate of the conductance autocor- 
relation to fit our experiments. We invoke the sim- 
plest non-equilibrium distribution for diffusing quasipar- 
ticles, that is the collisionless limit—, by lumping the 
relaxation processes in the damping parameter of the 
Cooperon/Diffuson propagators tc/d- The dependence 
on the applied voltage is introduced by tuning t v . We 
obtain oscillations in the negative tail of the autocorrela- 
tion, K g (V, 0, AV) (see Fig. (fit)]) ) and our scaling proce- 
dure fulfills the relation Dt v = L 2 . Extra contributions 
that are specific of the non equilibrium theory and are 
known to be responsible first for a linear increase of the 
autocorrelation function with V and subsequently for its 
power-law decay are not included here. 



VII. CONCLUSIONS 

We have reported about transport measurements of 
high quality biepitaxial Grain Boundary YBCO Joseph- 
son Junction at temperatures below 1 Kelvin, performed 
over a time period of about 18 months. A global 
view on the data offers a consistent picture, pointing 
to transport across a single SNS (superconductor nor- 
mal superconductor)-like diffusive conduction channel 
of mesoscopic size L < 0.1 fim. We have mostly ex- 
plored the magnetoconductance fluctuations in the volt- 
age range eV >> E c >> fc^T, where the Thouless energy 
E c ~ 1 meV . The Thouless energy, 2 — 3 orders of mag- 
nitude larger than the one usually experienced in nor- 
mal mesoscopic or low-T c superconducting samples, de- 
termines qualitatively the quantum coherent diffusion in 
the channel. We believe that the oscillations in the resis- 
tance that can be seen in Fig. ([3]) can be due to quantum 
diffusion. 

The mesoscopic correlations are found to be quite ro- 
bust in our GB narrow channel, even at large voltages. 
This could not occur if the lifetime of the carriers were 
strongly cut by non-equilibrium relaxation. We con- 
clude that mesoscopic effects deeply involve supercon- 
ducting electron-electron correlations, which persist at 
larger voltages. Transport features due supercurrents 
and quasiparticles at finite voltages cannot be disentan- 
gled in the pattern of the conductance, nor in its variance. 
This consideration has led us to approach the problem 
with a non-equilibrium model calculation for generic co- 
herent transport, which highlights the role of the phase 
breaking time t v , without including superconducting cor- 
relation explicitly. Fig l 101 shows the comparison between 
our model results and the autocorrelation experimental 
data, which is encouraging. The remarkably long lifetime 
of the carriers, which we find, appears to be a generic 
property in high-T c YBCO junctions as proved by optical 
measurements^ 2 - and Macroscopic Quantum Tunneling 13 . 
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